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1 Electromagnetic theory 

1.1 Photons and light 

Since the work of J.C. Maxwell in the late eighteen hundreds it is evident that light 
is most certainly electromagnetic in nature. In classical electrodynamics this leads 
to the picture of continuous energy transfer by means of electromagnetic waves. The 
modern quantum electrodynamic picture describes the transport of energy in terms of 
photons which are localized quanta of energy. This particle- wave duality will be used 
here to advantage. For the most part of what follows classical electrodynamics suffices 
to describe the reflection and refraction of light at surfaces. Ofcourse. quantum 
electrodynamics needs to be invoked when interaction of light and matter is treated 
(in section 1.3) but also when the origin of polarization is discussed (in section 1.5). 

1.2 Maxwell’s equations 

The earliest experiments on electromagnetic phenomena, that resulted in for instance 
Coulomb’s law for the force between two electrical charges and Biot-Savart’s law 
for the force between current carying conductors, are now captured by Maxwell’s 
equations supplemented with Lorentz’ expression for the force on a moving charge 
and with an expression for charge conservation. 

In Maxwell’s theory act an electric field E and a magnetic field B. Both are 
vector fields so that magnitude and direction of the electric field and the magnetic 
field may vary as a function of time and of position in space. Maxwell’s equations in 
fact summarize the conclusions of a great many of experiments on electromagnetic 
phenomena. It is exactly at this level where classical electrodynamics is empirical. 

1.2.1 Faraday’s induction law 

Faraday discovered that a time- varying magnetic flux, passing through a closed con- 
ducting loop, results in the generation of a current through that loop. When the loop 
is open, one measures the electromotive force that drives this current. Mathematically 
more sophisticated this is written as 

E -t dr = —dt j B ■ n d 2 r (1.1a) 

dA A 

where A is an open area bounded by the contour dA. Here and in the following 
t denotes the tangential unit vector associated with the line element dr and n the 
normal unit vector associated with the surface element d 2 r. The consequence of 
this equation is. that a time- varying magnetic field will always have an electric field 
associated with it. 

1.2.2 Gauss’ law, electric 

Another important observation is named after the German mathematician K.F. Gauss 
and relates the flux of the electric field through a closed surface 8V to the total charge 
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in the enclosed volume V. For a space charge distribution p{r ) Gauss’ law is written 
as 


/ 


E • n Sr 



(1.1b) 


The constant e is known as the electric permittivity and in free space its value is 
8.8542x 10 -12 C 2 N -1 m -2 in SI units 1 . Other media will be treated in section 1.3. 


1.2.3 Gauss’ law, magnetic 

There is no known magnetic counterpart to the electric charge. Magnetic fields are 
induced by electric currents, either macroscopic or microscopic as in ferromagnets, 
and hence the magnetic analog of eq.(l.lb) is 

(j> B -n d 2 r = 0. (1.1c) 

av 

1.2.4 Ampere’s circuital law 

As stated before, magnetic fields are generated by electric currents, that is currents 
of electric charges. This is the great finding of A.M. Ampere. Where Ampere had 
in mind currents as brought about in conductors by an electric field, Maxwell hy- 
pothesized the existence of displacement currents as occur between the plates of a 
capacitor while its charge is varied. Maxwell’s restatement of Ampere’s circuital law 
is therefore 

/B. t dr = fi j (j + ed t E) • n d 2 r (l.ld) 

8A A 

with j the current density. A is again an open area bounded by the contour dA. 
The quantity p, is called the permeability. For free space its value in SI units is 
47rxlO~ 7 Ns 2 C~ 2 . Other media will be treated in section 1.3. 

Just like eq.(l.la) the above equation has the important consequence, that a 
time-varying electric field always has a magnetic field associated with it. The two 
equations together lead to the wave equations that are discussed below. 

1.2.5 Maxwell’s equations in local form 

The set of equations (1.1a. .. l.ld) can be cast in a local (differential) form which 
makes them more useful for deriving the wave aspects of the electromagnetic field. 
This transition can readily be accomplished by making use of two theorems from 
vector calculus, namely Gauss’ divergence theorem 

f F-nd 2 r = Jv F d 3 r (1.2) 

dV V 


X A useful approximation is e 1/(367t10 9 ). 
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and Stokes’ theorem 

<j) F -t dr = j (V x F) ■ n d 2 r. (1.3) 

8A A 

with F an arbitrary vector field. The result is 

V x E = -d t B (1.4a) 

V • E = p/e (1.4b) 

V B = 0 (1.4c) 

VxB = p(j + edtE) (1.4d) 

In free space, where there are no charges ( p = 0) and where there are no currents (j = 
0), the combination of eq.(1.4a) and eq.(1.4d) immediately leads to two extremely 
concise vector equations 

V 2 E = epd/E (1.5a) 

and 

V 2 B = epd/B. (1.5b) 


These equations are the wave equations and they predict solutions in the form of 
electromagnetic waves with a velocity c = 1 / ^fep. It was exactly this prediction, 
which was experimentally verified by A.H.L. Fizeau. that made Maxwell’s formulation 
of classical electrodynamics one of the great intellectual triumphs of all times. 


1.2.6 General solutions 


The general solution of Maxwell’s equations in the presence of charges and electric 
currents can easily be formulated in terms of an electric potential <p(r, t ) and a mag- 
netic vector potential A(r,t ) such that 2 


and 


E = — V<£> — d t A. 


B = V xA 


Air ' t] = / 47r^|r — r’2 ( r ' T ^ ~ 



(1.6a) 

as 

(1.6b) 

d 3 r f 

(1.7a) 

d 3 r' . 

(1.7b) 


The time argument in the charge and current contains a term |r— r'|/c which accounts 
for the retardation: a disturbance at some location r' takes some time, even as it 
travels at the speed of light, to reach the position r. When retardation effects are 
not important this term can be neglected and Maxwell’s equations are treated in the 
static limit. 


2 There is some freedom in the choice of potentials but that is not of relevance here. 
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1.2.7 Multipole expansion 

The simplest electrostatic problem to consider is the case of a single charge located 
in the origin 

pi(r) = qS(r). (1.8) 


The electrostatic potential for this case is the well-known Coulomb potential 


Ti(r) 


Q 

Alter 


(1.9) 


The next simplest case is the physical dipole that consists of two opposite charges at 
a given seperation 


P2(r) 


< 1 (<*(»*- -Hr + ^)j 


with the electrostatic potential 


T2 (r) 


q 

Arte 



1 


\ 



( 1 . 10 ) 


( 1 . 11 ) 


The transition to a mathematical or point dipole occurs when taking the limit d — > 0 
while keeping the dipole moment p = qd constant: 


cp 2 (r) ~ 


1 p ■ r 

Aire r 3 


for d — y 0 while p 


qd constant. 


( 1 . 12 ) 


For what follows it is important to note, that the above potential can also be written 
as the derivative of the electrostatic propagator 


P2 (r) =p-V 


1 

47rer 


(1.13) 


Let us now consider an arbitrary charge distribution p(r ) inside a volume V 
around the origin and suppose that we wish to know the electrostatic potential in 
some far-away point r. By expanding the propagator in the electrostatic potential 



1 

Aire\r — r' 


p(r')dr'. 


(1.14) 


in a Taylor series. 

i 1 = - + r' • V- + ■p-.r'r' : VV- + . . . , 

T* — V\ T T 21 T 

the multipole expansion of the charge distribution is obtained 3 


(1.15) 


p{r) 


-\ q - 

Aite lr 


+ p- V- +Q : VV- + . . .} 

r r ) 


(1.16) 


3 Note that values of the multipole coefficients depend on the choice of the coordinate system 
except the lowest order non-zero coefficient. 
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with the total charge 

q = j p(r')dr' , (1.17a) 

v 

the total dipole moment 

p = j r'p(r')dr', (1.17b) 

v 

the total quadrupole moment 

Q = J r'r' p(r')dr' . (1.17c) 

v 

etc. Each term in this multipole expansion is a factor r~ 2 weaker than its predecessor. 
The conclusion of this analysis is therefore, that the lowest order non-zero multipole 
moment gives the dominant contribution to the electrostatic potential. 

1.3 Light and matter 

Although molecules are neutral a great many of them have a non-symmetrical charge 
distribution which gives rise to a permanent dipole moment. The water molecule for 
instance has a center of gravity of the positive charge that is more on the hydrogen 
side of the molecule whereas the center of gravity of the negative charge is more on 
the oxygen side of the molecule. It therefore carries a permanent dipole moment of 
6.17xlO -30 Cm. On the other hand the oxygen molecule is nicely symmetric and 
carries no permanent dipole 4 . In the absence of an electric field thermal agitation 
keeps the molecule randomly oriented so that there is a vanishing (time-) averaged 
dipole moment. With a (not too large) applied electric field E there is a tendency to 
line up along the field in the configuration of lowest energy leading to a nonvanishing 
average dipole moment 

<P> - 3 ( 1 - 18 ) 

with p 0 the permanent dipole moment. &#=1.38x 10 -23 J/K the Boltzmann constant 
and T absolute temperature, see exercise 1.1. 

The charges in a molecule are mobile and hence, when an electric field is applied, 
the electronic and even the atomic structure of the molecule is distorted and this, in 
turn, gives rise to an induced dipole moment 

p = aE + (3E 2 + .... (1.19) 

For not too large electric fields the induced dipole moment is linear in the electric 
field. The proportionality constant a is called the polarizability of the molecule and 
it is mainly of electronic origin. For higher electric fields the induced dipole moment 
is no longer linear in the applied field and this can be accounted for by for instance 
the hyperpolarizability (3 (see eq.(1.19)). 

4 Dipole moments of molecules, bonds, and molecular groups are sometimes expressed in Debije 
units: ID — 3.336 x 10 - 30 Cm. 
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1.3.1 Polarization 


Consider a capacitor consisting of two flat plates that have charges per unit area of +<r 
and —a respectively. In the absence of matter between the plates the magnitude of 
the electric field between the plates is equal to cr/eo where the dielectric permittivity 
of vacuum is denoted by £o When a (non-conducting) material is inserted between 
the plates Faraday already noticed that the magnitude of the electric field drops. For 
this situation one still writes for the magnitude of the electric field cr/e but now with 
a dielectric permittivity e > eo- The reason for the drop in the electric field is that 
the electric field effectively induces dipoles in the dielectric, either by aligning the 
molecular permanent dipoles or by inducing molecular dipole moments. For not too 
large electric fields the induced polarization is linear in the electric field 

P = e oX E, (1.20) 


where % is called the susceptibility of the dielectric. The magnitude of the electric 
field between the plates of the capacitor can still be related to both the surface charge 
and the polarization by 


E = 


a-P 

eo 


( 1 . 21 ) 


from which one derives e = eo(l + X ). Note that the polarization acts as a surface 
charge, see exercise 1.2. The quantity e r = e/eo = 1 + x is called the relative dielectric 
permittivity of the dielectric. 

The relation between the dielectric permittivity and the molecular polarizability 
of a dielectric is given by the Clausius-Mossotti relation 


( 1 . 22 ) 



e + 2e 0 


with p the molecular density of the dielectric, see exercise 1.3. This relation is derived 
under the assumption that the molecules in the dielectric are identical and uncorre- 
lated. It is trivially extended to the case of more than one species of molecules. Cor- 
relations are usually taken into account by means of a virial expansion. In particular 
for fluids (gases and simple liquids) the Clausius-Mossotti relation fits experimental 
data quite satisfactory. 


1.3.2 Dispersion 


Up till now only static electric fields were considered. For time- varying fields of low 
frequency the above relations still hold; this is the quasi-static regime. For higher 
frequencies the molecular susceptibility is frequency-dependent and this phenomenon 
is called dispersion. There are two distinct mechanisms that give rise to dispersion 
and both can be captured by the same formula 5 : 


a(co) = 


a 0 cjQ 


U!q + icncJoT — a; 2 


(1.23) 


5 In electromagnetism it is customary to use Fourier-transformed quantities. Here uj is to be 
regarded as the frequency of the applied electric field. 
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with cxq the low-frequency (static) limit of the polarizability. 

A permanent dipole rotates to follow the instantaneous direction of the electric 
field. Clearly, the rotational friction plays a dominant role but also the moment 
of inertia I of the molecule comes into play. By substituting q ; 0 = . 7 = 

£ r / y/2ksTI, and = 2ksT / 1 the response function eq.(1.23) accounts for this 
phenomenon. In many situations cn/cuo 1 and eq.(1.23) can be approximated by 
the well-known Debije relaxation formula 

(1.24) 

1 + IUT 


with the rotation correlation time r = 7 /cno = ^r/^ksT. The polarizability is a 
complex number. The real part is in phase with the incident field and the imaginary 
part is out of phase with the incident field. The imaginary part reflects the energy 
absorption: the rotational friction of the molecule causes electromagnetic energy to 
be transferred to the surrounding molecules. 

Each charge in a molecule can be considered as harmonically bound to its rest 
position. The restoring force is then proportional to cOq and the damping of the 
motion is measured by 7 . In principle the calculation of ao, uio, and 7 is a quantum 
mechanical problem but the thus desribed Sellmeier damped-oscillator model provides 
for an adequate picture, see exercise 1.4. 

In general there are more than one contribution to the dispersion and one writes 
for the polarizability 


k 


ul + iuJUklk ~ u 2 ‘ 


( 1 . 25 ) 


Typically there are three dispersion regimes: For low frequencies (static limit) the 
orientational polarization, the polarization arising from the permanent dipole mo- 
ments. is the most important. Since the rotation correlation time of a molecule is 
typically of the order of 1 ps this contribution vanishes in the microwave regime ( 10 9 
. . . 10 11 Hz). At the infrared frequencies ( 10 11 . . . 10 14 Hz) the distortion polarization 
vanishes. This polarization arises from distortions in the positions of the nuclei due 
to the electric field. At even higher frequencies only the electrons can still participate. 
At each absorption frequency, typically in the visible or ultraviolet regime, another 
electronic term in the polarizability vanishes. A typical graph of the polarizability as 
a function of frequency is given in fig. 1 . 1 . 


1.3.3 Magnetization 

For the magnetic properties of matter a presentation similar to the one in section 1.3.1 
can be given. In practice, however, it is hardly ever necessary to deal with magnetic 
properties of matter which is reflected by the fact that the relative permeability for 
most materials is equal to one. 


1.3.4 Refractive index 

As stated in section 1.2.5 electromagnetic waves, like light, have a velocity v = 1 / y/ejJ. 
This velocity is highest in vacuum, where it is approximately equal to 3xl0 8 m/s. 
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log (©) 

Fig. 1.1. Typical graph of polarizability versus frequency. 


In any other material this velocity is smaller and the ratio of the speed of an electro- 
magnetic wave in vacuum to that in matter is the absolute index of refraction of the 
material and it is given by 


n = 


efj, 


(1.26) 


Since for most materials the relative permeability is equal to 1 one has from eq.(1.26) 
for the relative dielectric permittivity e r = n 2 . Substituting this in the Clausius- 
Mossotti relation (1.22) with the polarizability given by eq.(1.25) yields the Lorentz- 
Lorenz relation: 


n 2 — 1 p 

— = — a 

n 2 + 2 3e 0 


(1.27) 


This relation shows that, for non-magnetic systems, the refractive index is frequency- 
dependent exactly in the same way as the dielectric permittivity. 


1.4 Maxwell’s equations in isotropic media 

In many materials, like most dielectrics, the induced polarization is indeed linear in 
the electric field so that eq.(1.20) holds. In a similar way one often has a magneti- 
zation M linear in the magnetic field. Both statements are ofcourse only true for 
a limited range of field strengths. For optical experiments this is almost always the 
case. It was for this reason that Maxwell used, besides the electric field E and the 
magnetic field B. the electric displacement D = e 0 E + P and the magnetic induction 

H = B/p 0 ~ M 6 . 


6 Sometimes one defines H' — B — poM such that II' and B have the same units. 
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The eqs.(1.4a. . . 1.4d) with the above mentioned substitutions then take a very 
simple form: 


V x E = -d t B 

(1.28a) 

V D = p' 

(1.28b) 

V B = 0 

(1.28c) 

VxH = j' + d t D 

(1.28d) 

which are to be used with the constitutive relations 

D = eE, 

(1.29a) 

and 

H = B/p. 

(1.29b) 


The charge density p' and current density j' are the charges and currents that are 
introduced into the medium. In free space the dielectric permittivity is eo and the 
magnetic permeability is p 0 . Substituting these for e and p in eqs. (1.28a. . . 1.28d) 
and in the constitutive equations (1.29a, 1.29b) yields Maxwell’s equations in local 
form (1.4a. . . 1.4d). 

1.5 Polarized light 

One of the simplest solutions of the wave equation (1.5a) is the monochromatic plane 
wave. The electric field of such a plane wave is given as 

E(r,t) = E 0 e exp {i(ut — k • r)} (1.30) 

where it is understood that, in order to obtain a physical field, the real part has to be 
taken 7 . E 0 is the amplitude of the electric field and e is a unit vector that determines 
the polarization of the light. Light is a transverse wave which means that the electric 
field varies in a plane perpendicular to the propagation vector k, i.e. e • k = 0. The 
magnitude of the propagation vector is the propagation number k = 2n/X with A the 
wavelength in the medium. It is related to the frequency co by the equation k = co/v 
where v is the wave velocity in the medium. It is customary to relate to the speed c 
of light and the wavelength Ao in vacuum: 

, 2irn con 

k = — = — 

Aq c 

with n the refractive index of the medium. 

The plane wave does not occur in nature but it is a convenient device for calcu- 
lations because any physical wave can be represented as a sum of plane waves. Two 

7 The magnetic field is always perpendicular to both the electric field and the propagation vector 
and. for nonmagnetic material, trivially related to the electric field. Therefore it is customary only 
to refer to the electric field. 


(1.31) 
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types of physical waves need particular attention. The one is the monochromatic 
spherical wave 

gi (ut—kr) 

E(r,t) = E 0 ^ (1.32) 

which is a good description for the light emitted by a point source. The other is the 
Gaussian beam which is a good approximation for a (single mode) laser beam 8 . A 
beam travelling in the ^-direction it is represented by 

l~2 e l Gt-kz) - (x 2 +y 2 ) /w 2 (z) 

E(r,t) = Ej — (1.33) 

V 7 r w 2 (z) 

with the position dependent spot size w(z) which, for 2 sufficiently large, is given by 
w(z ) ~ (1-34) 

TTW 0 

The waist w 0 is a design parameter of the laser. 

1.5.1 The nature of polarized light 

The plane of vibration, that is the plane in which the electric field varies, is per- 
pendicular to the propagation vector. An arbitrary coherent wave can therefore 
always be decomposed in two orthogonally vibrating component waves 9 . Assuming 
the propagation vector along the positive z-axis of a Cartesian coordinate system, 
these component waves can be represented as 

E x = 

(1.35) 

E y = E 0 ^ y e y e^- kz+ ^ 

where e x and e y are the unit vectors in the positive ^-direction and ^-direction 
respectively. E 0jX and Eq <y are the amplitudes of the electric field components and p 
is the phase difference. Only phase differences between — tt and 7r are relevant. 

Each of the components represents a linearly polarized wave. For p = 0 the two 
components are in phase and the resultant of the two components is also linearly 
polarized because the amplitude Eo iX e x + Eo ty e y is fixed. One could equally well have 
chosen the x-axis along the resultant and then the ^-component would have zero 
amplitude. For p = n/2 and equal amplitudes E 0jX = E 0ty = E 0 the resultant wave 
is left circularly polarized : the (real part of) the resultant field is written as 

E = Eq { e x cos(o )t — kz) + e y sin(aj£ — kz )} (1.36) 

and rotates counterclockwise in the plane of vibration, thus describing a circle. For 
p = — 7t/ 2 the wave is right circularly polarized. For arbitrary p the light is said to 
be elliptically polarized, see fig. 1.2. 

8 see A.E. Siegman. Lasers, Oxford University Press. Oxford, chapter 17. 

9 Coherent light consists of components that have a fixed phase difference in time. Natural light 
is incoherent and in practice light is only partially coherent. For most practical considerations laser 
light can be considered to be coherent. 
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Fig. 1.2. Various polarization states corresponding to specific values of cp (in degrees). 


The linear polarized wave seems to be the most straightforward from a mathe- 
matical point of view. Circular polarized light on the other hand seems to be a rather 
peculiar notion that can best be considered as the resultant of two linear polarized 
waves with n/2 phase difference. From a physical point of view, however, the opposite 
is the case! The polarization of light originates from the angular momentum of the 
photons: because of the quantized nature of angular momentum (a photon is a spin 1 
particle) photons can be considered as wave-packets which are either right-circular 
or left-circular polarized. Now linear polarized light appears to be the artifact that 
consists of a photon current, half of the photons in the left-circular state and the 
other half in the right-circular state! 


1.5.2 Jones notation 

As discussed before the plane of vibration is determined by two unit vectors. For 
polarized light it therefore suffices to discuss the projection of the electric field in the 
plane of vibration. This reduces the number of relevant components of the vector 
to 2 and such vectors are called Jones vectors after the american physicist R. Clark 
Jones who invented the notation in 1941 10 . 

To give an example for light that travels along the positive z-axis. the electric 
field vector is written as 

/ E 0 . x e^ 

E = E 0jy e i{py e i{ut ~ kz) (1.37) 

V o 

and the Jones vector is 


E 0 , x e llp 
E 0 , y e iVy 


(1.38) 


Each optical component has its own Jones matrix T) : and when light with Jones vector 
Ei passes through a sequence of these components the outcoming light has the Jones 
vector E 0 given by 
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where standard matrix multiplication is implied, see exercise 1.5. 
A horizontal linear polarizer has Jones matrix 

T = f 1 0 
Io ( 0 0 

and. using the matrix for rotation through an angle a 


(1.40) 


R 


a 


cos a sm a 
— sin a cos a 


(1.41) 


the Jones matrix for a linear polarizer of which the axis makes an angle a with the 
horizontal is calculated as 


T a — R- a T 0 R a 


cos 2 a 
cos a sin a 


cos a sin a 
sin 2 a 


(1.42) 


Another useful optical component is a retarder that effects a phase difference S be- 
tween its fast and slow axis 


Ts,o — 


e i5 / 2 0 

0 e - i<5 / 2 


(1.43) 


Retarders can be used to change for instance linear polarized light into circular polar- 
ized light and back. Although it is not often done it is possible to include polarization- 
dependent absorption in the Jones notation. 


Exercises 


1.1 The energy of a dipole moment p 0 in an electric field E is equal to«= — p 0 ■ E. Use 
the Boltzmann weight exp(— 0u) to derive 


ip) - Po 


(3poE smh.f3poE — cosh/3po-U ^ 

hj 

f3poE cosh {3poE 


with P = 1/kpT (Choose the ^-axis of the coordinate system along E). Subsequently 
show that for small relative energies fipoE this leads to eq.(1.18). 

1.2 Consider a parallel-plate capacitor with surface charge a. Calling the spacing between 
the plates d the voltage across the capacitor is given by the integral of the electric 
field over the distance between the plates. This yields 


V = 


ad 

eo 


and hence the capacitance is 


_ Q _ 6 od. 

~V~ ~d~' 
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where the plate area is taken to be A. 

(1) show that when a conductor of thickness b < d is inserted between the plates, the 
induced surface charge in the conductor leads to an increased capacitance 

_ e 0 A 

d{ 1 - b/d ) 

(2) when a dielectric is inserted the capacitance increases also due to an induced 
surface charge. This polarization charge arises from a net displacement 5 of the 
positive molecular charges with respect to the negative charges, see fig. 1.3. Calling 







± 

± 

± 

± 

± 

i 

L — 







i 

_ J 


Fig. 1.3. A dielectric slab in a uniform field. 

p the positive (or negative) charge density and q the elementary charge, the induced 
dipole moment is P — pqS. Argue, using fig. 1.3, that the polarization charge is equal 
to P and derive 

_ ep(l + x)A 
d 

using eq.(1.20). 

1.3 When an isotropic dielectric is placed in a homogeneous electric field E the molecules 
in the dielectric experience not the field E but a field that is modified by the presence 
of the other molecules. Let the polarization in the dielectric be P. Imagine we remove 
one of the molecules from the dielectric. What is left is a hole and for the electric 
field in the hole one derives 11 

^hole — E + - • 

3e 0 

The induced dipole in the molecule is then proportional to this field with polarizability 
a. For a density of p molecules per unit volume one has 

p= ^( e+ ^)' 

Use this to derive eq.(1.22). 


n see R.P, Feynman et al, Lectures on Physics II. sec. 11. 4, Addison- Wesley, Reading Mass., 1964, 
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1.4 Consider an electron, mass m and charge q. harmonically bound (force constant ko) 
to a nucleus (at rest) in an oscillating electric field E(t) — E$e lujt . Newton’s law gives 
the equation of motion of the electron in terms of the displacement x of the electron 
from the nucleus 

md 2 x — —kox + qE(t) 

The polarization of this simple system is P(t ) = qx(t). 

(1) Show that for frequencies uj ^ uq 

P(t) = ^z EocT*, 

where ujq — ko/m and ooccq = q/m. 

(2) The response at the resonance frequency uo is undetermined. Include a friction 
force —^dtx to overcome this problem and derive eq.(1.23). 

(3) Show that in the neighbourhood of the resonance frequency 

02 = qEo Q 

1 2mujQ (to — coo) 2 + Q 2 

where Q — £/2m. The plot of d 2 x versus frequency is called a Lorentzian distribution. 
It is commonly used to represent the shape of the lines in an absorption spectrum. 

1.5 Unpolarized light passes through a sequence of two linear polarizers of which the axes 
make an angle a with one another. Derive an expression for the transmitted light 
intensity as a function of the angle a: this is Malus* law. 
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2 Reflection at interfaces 

2.1 Geometrical considerations 



Fig. 2.1. Reflection and transmission of p-polarized light at a planar interface 
between two semi-infinite homogeneous optically isotropic media. 

Consider, with reference to fig. 2.1. the oblique reflection and transmission of a 
plane wave at the planar interface between two semi-infinite homogeneous optically 
isotropic media 1 and 2 with (complex) indices of refraction ni and n 2 , respectively. 
It is assumed that the change in index of refraction is abrupt across the interface, i.e. 
is a step function of the coordinate perpendicular to the interface. The wave incident 
in medium 1 gives rise to a reflected wave in the same medium and a transmitted. 
or more appropriately refracted wave, in medium 2. The angle of incidence 0*. the 
angle of refraction 6 t , and the angle of reflection 6 r are measured from the direction 
of the normal to the interface. The total fields inside media 1 and 2 obey Maxwell’s 
equations (cf. eqs. (1.28a . . . 1.28d)) and the boundary conditions at the interface. 
For the boundary conditions to be satisfied 

• the directions of propagation of the incident, reflected, and transmitted waves 
must all lie in the plane of incidence determined by the normal of the interface 
and the direction of propagation of the incident beam. 

• the angle of incidence must equal the angle of reflection, i.e. 

&i = 0 r , (2.1) 


and 


the angles of incidence and refraction must be related by Snell’s law: 
ni sindj = n 2 sin 6 t . 


( 2 . 2 ) 
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If media 1 and 2 are transparent, so that n i and n 2 are both real numbers, the angle 6 t 
is real. However, when either one or both media are absorbing, the angle 6 t becomes, 
in general, more complex and eq. (2.2) holds formally but the physical picture is more 
complicated. 

The law of reflection (2.1) is much more fundamental than it appears from this 
introduction. It can be derived easily by considering the interference of wavelets 
originating from the scatterers in the medium 1 . Consider as an example a homoge- 
neous distribution of scatterers in a plane, see fig. 2.2. A plane wave is incident on 



Fig. 2.2. Assembly of scatterers in a plane obeying the law of reflection. 

this plane with angle In all but one direction the spherical waves emerging from 
the scatterers interfere destructively. The one direction where there is constructive 
interference is in the plane of incidence and makes an angle 6i with the normal to the 
interface (see exercise 2.1). 

2.2 The Fresnel equations 

For a given amplitude and polarization of the incident wave, the amplitude and polar- 
ization of the reflected and transmitted waves can be determined from the continuity 
of the tangential components of the electric and magnetic field across the interface. 
It can be shown that when the incident wave is linearly polarized with the electric 
field parallel to the plane of incidence ( p-waves ), the reflected and transmitted waves 
are also linearly polarized with their electric field vectors parallel to the plane of in- 
cidence. Likewise, when the incident wave is linearly polarized perpendicular to the 
plane of incidence ( s-waves 2 ), the reflected and transmitted waves are also linearly 
polarized perpendicular to the same plane. 


1 For crystals the phenomenon is known as Bragg reflection. 

2 The letter s originates from the German word senkrecht for perpendicular. 
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In terms of Jones notation (cf. section 1.5.2) the transfer matrix for reflection or 
transmission has two eigenvectors corresponding to linear polarized light with electric 
field parallel (p) and perpendicular (s) to the plane of incidence. This is. however, 
only true if both media are optically isotropic as is assumed in what follows. The 
Jones matrix for reflection (transmission) is. in general, not energy conservative since 
part of the energy is transmitted (reflected). Let (Ei p ,E is ). ( E rp ,E rs ), and ( E tp ,E ts ) 
represent the Jones vectors for incident, reflected and transmitted light respectively, 
then the Jones matrix for reflection is 

*=(or.) < 2 - 3a) 


and for transmission it is 


t p 0 
0 t s 


(2.3b) 


The reflection and transmission amplitude ratios, due to Augustin Jean Fresnel (1788 
- 1827), can be derived by matching the tangential electric and magnetic field across 


the interface. They are 3 


n 2 cos 6i — ni cos 9 t 

r p — 

77.2 cos 6i + ni cos 9 t 

(2.4a) 

ni cos 6i — n 2 cos 9 t 

r s — 

77i COS + 772 COS 9 t 

(2.4b) 

2 Tii COS Oj 

t — 

p n 2 cos 9i + 77 1 cos 6 t 

(2.4c) 

2ni cos 9j 

t s = 

77 1 COS 9i + 772 COS 6>t 

(2.4d) 


The incidence angle 9i and the angle of transmission 9 t are related by Snell’s law 

( 2 . 2 ). 

The reflective properties of interfaces are described either by the reflectivities 
R p = | r p | 2 and R s = \r s \ 2 for p- waves and s- waves respectively, or by the ellipsometric 
angles ij} and A determining the (complex) ratio 

— = tan ipe' A (2.5) 

r s 

In figure 2.3 two examples of reflectivities are given, one for the air/glass and the other 
for the glass/air interface. The marked difference between the two is the angle of total 
reflectance for the glass/air interface which is absent for the air/glass interface. The 
existence of such an angle can be deduced from Snell’s law (2.2): upon transition 
from a dense into a less dense medium the refractive index jumps to a smaller value 
and hence the angle of propagation with the normal increases. Increasing the angle 
of incidence, the angle of refraction in the less dense medium reaches 90 degrees 
sooner than does the angle of incidence. Beyond this angle of total reflection no 
light is refracted and all is reflected. In actual fact, some electromagnetic field is still 
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present near the interface in the other medium, this is called the evanescent wave. At 
normal incidence (0j = 0) there is no distinction between p-waves and s-waves and 
consequently there can be no difference between the normal reflectivities. The p-wave 
reflectivity passes through a minimum. R p = 0, at the Brewster angle or polarizing 
angle. Inspection of eq. (2.4a), with the help of Snell’s law (2.2), shows that this 
angle is given by 

6b = arctan ^ • (2-6) 

At this angle the incidence wave is totally refracted into the second medium. Also, 
the p-wave amplitude reflection coefficient (2.4a) changes sign at this angle which 
causes the ellipsometric angle A to jump 180 degrees, see fig. 2.4. 

2.3 Layered structures 

The simplest layered structure where a film of refractive index no coats a substrate of 
refractive index n\ is shown in fig. 2.5. The medium above the film has refractive index 
U 2 and furthermore the three media are considered to be isotropic 4 . The procedure to 



Fig. 2.5. Oblique reflection and transmission of a plane wave by a layered 
system with parallel-plane boundaries. 

be followed to obtain the effective amplitude reflection and transmission coefficients 
is due to Airy 5 and will be left as an exercise (2.2). The result can, both for p-waves 


and s-waves, be written as 


_ no T r Q 2 e 2iip 

1 + r w ro2e 2i,p 

(2.7a) 

t _ t 10 t 02 e l * 

1 + r w ro 2 e 2 ^ 

(2.7b) 


where and t ,. ; are the amplitude reflection and transmission coefficients respectively 
for the i to j interface, see eqs. (2.4a. . . 2.4d). The phase angle ip arises from the 
optical path difference between the directly reflected light and the light that has 

3 The directions of p and s are chosen such as to form a right-handed Cartesian coordinate system 
with the direction of propagation. 

4 The more complex cases can be found in R.M.A. Assam and N.M. Bashara. Ellipsometry and 
Polarized Light, sec. 4. 7, North Holland, Amsterdam, 1977. 

5 G.B. Airy, Phil. Mag. 2 (1833) 20. 
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undergone one transmission into the film, one reflection and then again a transmission 
out of the film: 

27 m 0 d 

ip = COS Uq 

A 

Note that Snell’s law still holds, i.e. 

Tii sin B\ = Tig sin Oq = ri 2 sin 0% • (2-9) 

Since the films are usually thin, i.e. d A, the eqs. (2.7a) and (2.7b) are single- valued 
in d. For thick films this is no longer the case. The reason for this is in eq. (2.8): the 
phase angles cp and cp + it result in identical reflection and transmission amplitude 
coefficients. 

2.4 Thin films 

The ambient/ film/ substrate system is often used as a model system in experimental 
situations. The layered structure of the previous section 2.3 is one of the crudest 
examples since the refractive index will in practice vary much less abrubtly then at 
the boundaries, see fig. 2.6. 


(2.8) 


£ 



Fig. 2.6. Example of a dielectric permittivity profile with limiting values ei 
and 62 far away from the interface. 

Several methods have been developed in order to deal with more realistic systems. 
One method consists of extending the layered system of section 2,3 to a multilayer 
system upon which a limiting procedure is performed. In the limit of infinitesimally 
thin layers a continuous refraction index profile is obtained 6 . Another method treats 
the film as a perturbation to the ambient /substrate interface. It has been developed 

6 see R.M.A. Assam and N.M. Bashara. Ellipsometry and Polarized Light sec. 4. 6. North Holland. 
Amsterdam. 1977. 
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by Lekner for stratified media 7 and independently by Bedeaux and Vlieger for general 
surfaces 8 . 

For thin films the second method is the more appropriate because interface rough- 
ness can also be taken into account. Therefore this method is discussed below al- 
though for pedagogical purposes the system will at first be assumed to be stratified, 
i.e. the dielectric permittivity varies only in the ^-direction perpendicular to the in- 
terface. At the end of this section surface roughness will be discussed briefly. Both 
ambient and substrate are supposed to be isotropic and all three media are non- 
magnetic. 

In the absence of a film, the dielectric permittivity profile is given by 9 

eo {z) = ei 6{—z) + e 2 d(z). (2.10) 


This profile gives rise to the Fresnel reflection and transmission amplitude coefficients 
discussed in section 2.2. An important aspect of the method is that not all of the 
details of the profile e(z) need to be taken into account. For sufficiently thin films a 
very limited number of moments of the dielectric permittivity profile suffices. These 
moments are defined as integrals of the excess of the dielectric permittivity with 
respect to the Fresnel profile (2.10). Some of them are given by 


3 


3 — 



k = 0,1,2,... (2.11) 


there are however more complicated moments. From the above definitions it is clear 
that the k-th moment is of k-th. order in the thickness of the film 10 (with respect 
to the wavelength of the light). Another aspect is that the location of the origin of 
the coordinate system manifestly modifies the values of the moments. Obviously, the 
measurable quantities such as the reflectivity are independent of the location of the 
coordinate system. These are expressed in invariant combinations of the moments. 
One first order invariant moment is for instance constructed as J7i = j+,i — f -i^ 23 -,i- 
For non-absorbing films the three lowest order invariant moments are: one of first 
order in film thickness 


Ji 



fii - e(z)}{e(z) - e 2 } d ~ 
e(z) 


( 2 . 12 ) 


and two of second order in film thickness 

j21 = (p 2 // ^ ~ 1 ^ ~ 2) W 

7 see J, Lekner. Theory of Reflection, chapters 3 and 8. Martinus Nijhoff Publ.. Dordrecht 1987. 
8 see e.g. D. Bedeaux. G.J.M. Koper. E.A. v.d. Zeeuw. J. Vlieger. and M. Wind. Physica A 207 
(1994) 285 and references therein. 

9 0(z) is the Heaviside function: it is 1 for z > 0 and it is 0 for z < 0, 

10 The thickness of the film is size of the region in which the refractive index profile e(z) differs 
from the Fresnel profile eo(^). 
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X { f (y) e 0 (2'-z)} <2 ' 13) 

J 22 = - (f)7/ {fW ~ £0( " ~ ~ C ° (2 ' ~ Z)) izdz'. (2.14) 

Absorbing films, however, have a complex dielectric permittivity and therefore also 
have two imaginary invariants of first order. Im(j_i) and Im(j +1 ). For simple profiles 
these invariants are easily calculated, see exercise 2.3. 

For the reflectivity for s- waves for an absorbing film, the lowest order contribution 
comes from a first order imaginary invariant moment: 


R s 


K S}0 ~ 4im(j_ 




(2.15) 


where for non-absorbing films the lowest order contribution comes from a second 
order invariant moment 


R. 


Rs, 0 


2 ei J22 


4nin 2 cos cos 0 2 
(n 2 cos 0i + rii cos 0 2 ) 4 


(2.16) 


For the reflectivity for p-waves for an absorbing film the lowest order contribution is 
also from a first order imaginary invariant moment: 


R n 


0 “T 


X 


{(^€2 sin 


(n 2 cos + 
2 0i Im(j +; i) — cos 


2 02 Im(j 


(2.17) 


where for non-absorbing films the lowest order contribution is of second order in film 
thickness 


Rn = 




(n 2 cos 0i 


~ 1 


X 


-2(J7 2j i - €1 sin 2 Q\J 2 ^){e 2 cos 2 0 1 - e x cos 2 0 2 ) + — sin 4 0 X 

e 2 J 


(2.18) 


The values R Sj0 = ks,o| 2 an< i R p , 0 = \r p ,o\ 2 are those found for the Fresnel profile 
eo(z). Contrary to the reflectivities, the ellipsometric ratio, cf. eq. (2.5), already has 
a contribution of first order both in the absorbing and in the non-absorbing case 11 


r JL 

r s 


r s , 0 


2 i ni cos 0i sin 2 0i 
r s,o ( n 2 cos 0i + n 1 cos 0 2 ) 2 


Ji. 


(2.19) 


The case of reflection by a thin film between like media where e x = e 2 , like soap films, 
requires special attention 12 . For the Fresnel profile the Brewster angle is defined by 
eq. (2.6). In ellipsometry it is common practice to define the Brewster angle as the 
angle for which the ellipsometric angle A (see eq. (2.5)) vanishes. The shift in the 

1:L This formula is often attributed to Drude, but was, in essence, first obtained by L. Lorenz in 
1860. 

12 see for instance J. Lekner .Theory of Reflection, chapter 3. Martinus Nijhoff Publ., Dordrecht 
1987. 
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position of this Brewster angle due to the film is of second order in the film thickness. 
The contribution to the ellipsometric ratio at the Brewster angle is (for the Fresnel 
profile it vanishes) 


(0b) 


\/ e 2 + e l 
2 1 €2 — Ci 


Ji- 


(2.20) 


For non-absorbing films, as discussed above, the contributions to the reflectivities 
and the ellipsometric ratio is small. The notable exeption is the p-wave reflectivity 
and the ellipsometric ratio near the Brewster angle. Therefore most experiments are 
carried out near this Brewster angle as will be discussed in section 2.5. 

The important result of the work of Bedeaux and Vlieger is that the functional 
form of the reflectivities (2.15, 2.16, 2.17, 2.18) and the ellipsometric ratio are pre- 
served when interface roughness is introduced. The interpretation of the invariant 
moments however is completely different. For instance for a film consisting of a 
low density 1 ? of small dielectric spheres (radius a. polarizability a) Haarmans et al. 


calculate 13 



Ji ^ 

444 ) 

(2.21) 

J21 ^ 

V A / 

(2.22) 

J22 ^ 

4 ( 4 - 

(2.23) 


where it is assumed that the dielectric permittivity of the spheres is close to the 
permittivities of both the medium and the ambient. More exact formulae can be 
obtained from the cited paper. 


2.5 Experimental techniques 

In this section three techniques that can be used to study the reflection from interfaces 
will be discussed. 


2.5.1 Null ellipsometry 

In its simplest form a null ellipsometer consists of two telescope arms (see fig. 2.7) 14 . 
One arm holds a light source such as a small laser, a polarizer, and a linear retarder. 
This arm directs the light to a reflecting surface, the sample. The angle of incidence 
can be chosen by means of a goniometer that rotates either the arm or the sample. 
The other arm, holding an analyzer and a photodetector, can be rotated around the 
sample in order to direct the reflected light from the sample into the photodetector. 

To analyze this setup we assign angles 9 p ,6 r , and 9 a to the orientation with respect 
to the plane of incidence of the fast axes of the polarizer, the linear retarder, and the 
analyzer respectively. After the polarizer the light, which is assumed to be of unity 

13 M. Haarmans and D. Bedeaux, Thin Solid Films 224. 1994, 117, Physica A 207 (1994) 340. 
14 see R.M.A. Assam and N.M. Bashara, Ellipsometry and Polarized Light, chapter 3. North Hol- 
land, Amsterdam, 1977. 
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Fig. 2.7. Null ellipsometer with on one arm the laser (L), the polarizer (P), 
and the linear retarder (R). In the middle the surface under study (S), and on 
the other arm the analyzer (A) and the photodetector (D). 


intensity, has Jones vector (cos8 p ,sm8 p ). The transfer matrix of the linear retarder 
is written as 


1 0 
0 p r 


(2.24) 


where p r is simply i for an ideal quarter wavelength retarder. Calling the amplitude 
reflection coefficients of the sample r s and r p , the electric field of the light impinging 
on the photomultiplier for unity incident intensity is given by 

( cos 2 6 a cos 6 a sin 8„\ ( r„ 0 \ 

E ° = (oo.fl.BinJ. «n’«. )(or.) 

I cos 2 6 r T p r sin 2 8 r (1 — p r ) cos8 r sm8 r \ ( cos0 p 
l (1 — Pr) cos 8 r sin 8 r sin 2 8 r + p r cos 2 8 r J l sin 8 P 



The condition(s) for which the detector signal vanishes gives the ellipsometric ratio 


r p tan0 r + p r tan(0 p — 8 r ) 

r s a 1 — p r tan 8 r tan(# p — 8 r ) 


(2.26) 


An advantage of the technique is that this result does not depend on the intensity 
of the incident light. Of course the accuracy with which the nulling angles can be 
obtained does depend on light intensity. Null ellipsometry can be used in cases where 
the film (or the substrate) is absorbing. As discussed in the previous section, the 
contribution of dielectric films to the ellipsometric ratio is small except near the 
Brewster angle. But even then the ellipsometric ratio is smaller than 10 -3 (it is 0 for 
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the Fresnel profile). Small errors in the nulling angles introduce too large an error in 
the ellipsometric ratio to use this technique succesfully for these films 15 . 


2.5.2 Polarization modulation ellipsometry 

There are two situations where nulling ellipsometry is not the technique of choice. 
The one situation is where the sensitivity of the technique is insufficient as discussed 
before and the other situation is where speed is needed in order to follow for example 
an adsorption process. In those situations polarization modulation ellipsometry is the 
better choice since it is fast, accurate and less sensitive than nulling ellipsometry to 
the alignment. The technique is used by Meunier 16 after the design of Jasperson and 
Schnatterly 17 and Beaglehole 18 . The optical arrangement is almost identical to the 
nulling ellipsometer, see fig. 2.7. It has no linear retarder and there are two additional 
elements, a polarization modulator operating at a high frequency of approximately 
50 kHz between the polarizer and the sample and one operating at a low frequency 
(of approximately 10 Hz) between the sample and the analyzer. The polarizer is 
oriented at 45° with respect to the plane of incidence so that the electric field has 
equal s and p amplitudes. The first modulator introduces a time dependent phase 
shift 5 = 5osino;t between the s and p components. The phase amplitude 5o can be 
controlled by the amplitude of the voltage applied to the modulator. After reflection 
at the surface the light passes through the second modulator that is operated at a 
low frequency. This modulator is calibrated such that it introduces a phase shift 
of either 0 or ir radians between the p and the s component at the low switching 
frequency. Then the light passes through the analyzer that is oriented at 45° with 
respect to the plane of incidence and impinges onto the photodetector. The action of 
both the second modulator and the analyzer is that of an analyzer switching between 
+45° and —45° with respect to the plane of incidence. Usually the photodetector is 
equipped with a pinhole to reduce the influence of stray light 19 . 

The analysis of this configuration starts by writing down the Jones vector of the 
electric field at the detector for unity incident intensity 


Ed — 


1 

2y/2 


1 ±1 
±1 1 


r p 0 
0 r s 



(2.27) 


and after some algebra one finds that the intensity EdE ^ at the detector is propor- 
tional to 


Id ~ H 2 ± 2|p| cos(A + <5) + 1 (2.28) 


15 There is however the notable exception of C. Bouhet. Ann Phys (Paris) 15 (1931) 3. who used 
the technique for surfactant films! 

16 J. Meunier. J Phys (Paris) Lett 46 (1985) L-1005, 

17 S.N. Jasperson and S.E. Schnatterly. Rev Sci Instr 40 (1969) 761. 

1S D. Beaglehole. Physica B 100 (1980) 163, 

19 The size of this pinhole is not always unimportant as it determines the collection angle for 
scatter contributions from the interface, see for instance J. Meunier. J Physique 48 (1987) 1819 and 
E.M. Blokhuis and D. Bedeaux. J Chem Phys 95 (1990) 13, 
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where \p\ is the absolute value of the ellipsometric ratio and A its phase. As 
5 = So smut the detector signal will contain besides a constant (0) component also 
components at frequency u, 2 u>, etc. The first few components are given by 

To ~ \p\ 2 ±2|p|J 0 (Ao)cosA + 1 

I w ~ ±4|p| Ji(5o) sincatsin A (2.29) 

7 2w ~ ±4|p| J 2 (Jo) cos 2o>tcos A. 

where J n is the Bessel function of order n. The procedure is now as follows: the low 
frequency switching of the signal makes it possible to find the phase amplitude Sq 
at which Jo(5o) = 0: then To remains constant while the modulator switches. The 
value of do at which this occurs is precisely known and therefore also the values of 
MSq) and J 2 (^o) are known constants. From the remaining signals 7 0 , 7 W , and 7 2al 
the ellipsometric ratio can be determined. 


2.5.3 Ellipsometry data analysis 

Both ellipsometry techniques provide the ellipsometric ratio. Without an underlying 
model this (complex) number is useless. In some instances one does have a model, 
such as the capillary wave model for fluid/fluid interfaces that will be discussed in 
the next chapter 3. More often there is no model at hand and one takes the default 
thin film model. The best situation is when both the ambient and the substrate are 
well characterized, i.e. the refractive indices are known and the uncoated interface 
is sufficiently smooth so that it does not act as a surface film itself. When these 
conditions are fulfilled, the ellipsometric ratio determines the (complex) first invariant 
moment ffi through eq. (2.19). Again, when no information is available on the 
dielectric permittivity profile one usually takes the simplest case of a homogeneous 
film of dielectric permittivity e with a thickness d. For such a film the first invariant 
is 

J, = (2.30) 

A 6 

The choice of the dielectric permittivity determines the thickness d that one will 
find: one finds an ellipsometric thickness. When a homogenous film of some material 
is adsorbed one may choose the dielectric permittivity of the bulk material. For 
discontinuous island films a method by Maxwell Garnett 20 is rather popular. Calling 
eo the dielectric permittivity of the film material and <p the volume fraction taken by 
the material rather than by the ambient in the film the effective dielectric permittivity 
e is found from 


e — e 2 _ e 0 - e 2 
e + 2e 2 €q + 2e 2 


(2.31) 


a formula that can easily be derived from the Clausius-Mossotti or Lorentz-Lorenz 
relation (1.27). A better approach to discontinuous films is the sphere model discussed 
at the end of section 2.4. 

20 see R.M.A. Assam and N.M. Bashara, Ellipsometry and Polarized Light, section 4.8, North 
Holland, Amsterdam. 1977. 
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An important application of ellipsometry is spectroscopy where the wavelength 
of the light is varied. With this technique, in particular when infrared light is used, 
it is possible to identify the adsorbed molecules at the interface. Even the type of 
binding to the substrate can sometimes be determined. 

2.5.4 Reflectometry around the Brewster angle 

The main drawback of ellipsometry is that it is practically impossible to obtain more 
information on the interfacial region than the first invariant moment. Even the de- 
termination of the ellipsometric ratio at a number of different incidence angles does 
not really improve this situation because usually the first invariant will dominate. 
This is not the case for p-wave reflectivity of dielectric films, see eq. (2.18). When 
p-wave reflectivities are collected near the Brewster angle, where the Fresnel contri- 
bution R p )0 is small, both second invariant moments and the first invariant moment 
(squared) can be determined 21 . 

The optical arrangement of a reflectometer is even simpler than the null ellip- 
someter, see fig. 2.7: there is no linear retarder 22 . The simplicity of the instrument 
allows for a variation of the incidence and detection angle. The technical difficulty is 
to keep the telescope arms in the plane of incidence of the reflecting interface. 

Exercises 

2.1 Argue that a random assembly of scatterers all lying in the same plane obey the law 
of reflection. 

2.2 Derive eqs. (2.7a, 2.7b) for the reflection and transmission amplitudes of a bilayer. 

2.3 Calculate the invariants Ji, J 21 , and J 22 for the homogenous film of dielectric per- 
mittivity e and thickness d. 

2.4 In a poor man’s ellipsometer 23 one measures R p / R s rather than the ellipsometric ratio 
r p /r s and hence one determines J7 2 (see eq.(2.19)). Demonstrate that the sensitivity 
of the method can be greatly improved by pre-adsorbing a layer of given ellipsometric 
thickness Jip. 


21 Note that the p-wave reflectivity is linear in j} , S 21 and Jii- 
22 see G.J.M. Koper and P. Schaaf. Europhys Lett 22 (1993) 543. 

23 see J.C. Dijt, M.A. Cohen Stuart. J.E. Hofman, and G.J. Fleer. Colloids and Surfaces 51 (1990) 
141 - 158. 
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3 Reflection at fluid/fluid interfaces 

3.1 The capillary wave model 

Consider the interface between two (immiscible) fluids. Below the interface the den- 
sity is pi and the dielectric permittivity is ei, and above the interface these are p 2 and 
e 2 respectively. A coordinate system is chosen such that the 2 -axis is perpendicular 
to the interface. The separation between the two fluids is a corrugated interface £(*) , 
where a? is a (2-dimensional) vector in the xt/- plane. Across this interface all physical 
properties change abruptly from their one value to the other. The interface does have 
a surface tension and a bending elasticity that try to keep the interface flat. However, 
due to thermal agitation it fluctuates constantly around its mean value (£) = 0 such 
that the interface appears to have a finite thickness (( 2 ). In the following section 
an expression for (( 2 ) will be derived and in section 3.1.2 the reflectivity of such an 
interface will be calculated. 

3.1.1 Statistical thermodynamics 

The statistical mechanical treatment consists of replacing the time-averages of the 
fluctuating interface by the associated ensemble averages over all possible corrugated 
interfaces weighted by their free energies. 

As a reference state the flat interface, Cl*) = 0 is chosen. The first contribution 
to the free energy of an arbitrary corrugated interface C(*) is the hydrostatic energy 1 

-Fhydr = ^ j A pg C 2 (*) d 2 x, (3.1) 

s 

where A p = pi — p 2 and g is the gravitational accelaration. The second contribution 
comes from the increase in surface area with respect to the reference state 

AA = /{Vl + |V C (*)| 2 -l} d 2 x~^j |V((*)| 2 d 2 x, (3.2) 

s s 

where it is assumed that the gradients in the corrugations are small 
free energy is therefore 

F surf = 7 AA = || |VC(*)| 2 d 2 x, 
s 

where 7 is the surface tension. The third contribution comes from the bending of the 
interface with respect to the reference state. For small corrugation curvature |V 2 C| 
it is given by 2 

F curv - \ J |V 2 C(*)| 2 d 2 x, (3.4) 

5 


. The increase in 

(3.3) 


1 The surface integrals extend over a macroscopic surface S 
2 see W, Helfrich. Z Naturforsch C 28 (1973) 693, 
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where /c is the bending elasticity modulus. In principle it is possible to include more 
contributions to the free energy of a corrugated interface. It is believed that the three 
mentioned contributions are sufficient to model real liquid/liquid interfaces. 

At this stage it is useful to introduce 2-dimensional Fourier-transformed quantities 
£(q) such that 

C(*) = f C(<?) exp (i q x) d 2 q. (3.5) 

Each value of q represents a capilary wave or mode, hence the name of the model. 
The total free energy of a corrugated interface £(aj) is. with eq.(3.5), written as 

F = 2tt 2 | {Apg + -fq 2 + ^}\aq)\ 2 d 2 q. (3.6) 

In order to find the mean squared thickness the thickness of all corrugated interfaces, 
weighted with their appropriate Boltzmann factors, have to be averaged. This is 
identical to averaging the thicknesses of all modes with the appropriate Boltzmann 
factors. Since each mode represents one degree of freedom the equipartition theorem 
states that the free energy associated with each such mode equals k B T / 2: 


2\ k B T 

(27r) 2 (Ap g + yq 2 + nq 4 ) ' 

The mean squared thickness of the interface is now calculated as 


(3. 

with the capillary wavelength L 2 = 7/A p g and the bending length L b = /c/7. The 
maximum wavenumber (/ lliax is determined by experimental conditions. In particular 
it cannot be larger than the inverse molecular spacing. When both the capillary 
wavelength L c and the bending length L b are not too small the resulting mean squared 
thickness is independent of g max . The first length to become too small is the bending 
length L b . Then the resulting mean squared thickness depends on the maximum 
wavenumber g max 3 . 

For too large corrugations and corrugation gradients the linear approximations 
that have been made to derive the free energy of a corrugated interface are no longer 
correct. For sufficiently small nonlinearities a mode-coupling theory has been devel- 
oped by J. Meunier 4 . 


<c 2 > = 



k B T 

4tt7 


log' 


1 + Qm&xF 


dq 


1 + 9i 


2 r 2 

max . 



( 3 . 7 ) 


3 see J. Meunier. J . Phys. Lett. (France) 46 (1985) 1005-1014, 

4 see J, Meunier. J. Phys. ( France ) 48 (1987) 1819-1831, 
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3.1.2 Reflectivity 

For quite some time it was believed that the finite thickness (£ 2 ) of the corrugated 
interface was the important contribution to its reflectivity. It is still instructive to 
repeat the argument as the exact theory is too complicated to render in these notes 5 . 
As a first step the averaged interfacial composition is calculated. For simplicity it is 
assumed that the two fluids are two phases of one and the same type, such as a liquid 
and a gas phase: 

p(z) = Pi + (to - Pi)(0(z - C(x))) (3.9) 

The average of the Heaviside function is calculated using the identity 

U 

0{u) = / S(u') du ' (3.10) 

— OC 

and by subsequently introducing the Fourier-integral representation of the delta- 
function 

S(u) = — f e iau da. (3.11) 

2tt J 

Substitution of these identities in eq.(3.9) yields 

p(z) =pi + (p 2 - Pi) J — Jda e lcnt (exp (-iaC(ai))) du. (3.12) 

— OC 

The average of the exponential is calculated with the following “trick” 

(exp(-iaC(*))) = explogjl - io(C) - |a 2 (C 2 ) + • ■ ■} 

= exp {-ia(C) - |a 2 ((C 2 ) - (C) 2 ) + ■••}■ (3.13) 

Since (() = 0 the leading contribution is due to the thickness (£ 2 ): 

(exp (-iaC(aj))) ~ exp j-y (C 2 >| (3.14) 

Finally one obtains from eq.(3.9) the result 


/x p2 T p2 . Pi p2 r 
P\ z ) = — 2 — + — 2 — erf 

with the error function defined as 



\d<d> 


(3.15) 


(3.16) 


Through the Lorentz-Lorenz equation (1.27) the averaged dielectric permittivity 
profile e(z) is related to the averaged density profile (3.15). An almost identical 


5 see F.P. Buff, R.A. Lovett, and F.H. Stillinger Jr. Phys. Rev. Lett. 15 (1963) 621. 
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result is obtained by using the (simpler) density-weighted average of the dielectric 
permittivity 


e(z) 


€l{Pl ~ P( Z )} + e 2 {p{z) - p 2 } 
Pi — P2 


Cl + €2 
2 + 



\/2<0 


( 3 . 17 ) 


At the Brewster angle the ellipsometric ratio is determined by the first invariant 
moment (2.12) which for this situation is given by 6 : 


Ji 


'' '-1 


27T e 1 — €2 


A 6i + €2 




1 — erf 


1 + ^ 


( 3 . 18 ) 


For the thickness (£ 2 ) the expression (3.8) can be used. This formula is not correct as 
it does not have the temperature dependence and the dependence on surface tension 
as have been found by experiment. The above given derivation is wrong in neglecting 
the effect of the slopes of the corrugated surface. These introduce phase differences 
between the p and the s waves that can only be properly accounted for by applying 
Maxwell’s equations to the interfacial regime. The result of such an analysis is 7 


Ji = 


3tt (d - c 2 y 


1 9<IC(9)| 2 } 


d 2 q. 


(3.19) 


A Ci + c 2 

The integral over q is easily evaluated in the limit of the capillary length L c — > oo 8 : 


/ 9<IC(9)| 2 > 


A 


k B T 


arctan(<? max L 6 ) 


27t^/tk 

For not too large a bending length L& the result is independent of q max : 


(3.20) 


37 r (ei - e 2 ) 2 k B T 

4A + e 2 ^/y/5 


(3.21) 


The thickness (£ 2 ) does play a role in the reflectivity of corrugated interfaces: the 
two second invariant moments (see eqs.(2.13,2.14)) are proportional to this thickness. 

6 see D. Beaglehole. Physica B 100 (1980) 163, 

7 see B.J.A. Zielinska. D. Bedeaux. and J. Vlieger. Physica A 107 (1981) 91. 

8 see E.M. Blokhins and D. Bedeaux. Physica A 164 (1990) 515, 
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3 REFLECTION AT FLUID/FLUID INTERFACES 


3.2 Application 

The direct application of the capillary wave model is hampered by the fact that 
realistic interfaces are not ideal, i.e. the dielectric permittivity does not step from 
one into the other value across the interface. From their general theory for rough 
surfaces Bedeaux and Vlieger find that if several surface inhomogeneities are present, 
the nett invariant moments are the sum of the separate contributions. So the above 
derived roughness term is at least to be added to a term due to the intrinsic density 
profile. 

For the interface between two fluid phases the classical density profile of the Van 
der Waals and Cahn-Hilliard theories can be used 9 

P( z ) = Pl ^ p2 + Pl 2 P2 tanh ^ . (3.22) 

This density profile is to be substituted in the expression (2.12) for the first invariant 
moment using the Lorentz-Lorenz relation (1.27). Far from the critical point the 
correlation length £ is small and the interface is almost like a step function. But as 
the critical point is approached, the interface thickens and becomes diffuse, until, at 
the critical point, it disappears completely and the two phases become identical. More 
complicated interfaces require a more detailed model for the dielectric permittivity 
profile. 

The apparent dependence of the ellipsometric ratio at the Brewster angle on 
bending elasticity modulus and surface tension as expressed by the first invariant 
moment (3.21) has been utilized by Meunier to determine the bending elasticity 
of surfactant monolayers at the planar oil/water interface 10 . By varying the salt 
concentration of the water phase samples with different surface tension and the same 
bending elasticity are prepared. The surface tension is determined independently. 
The obtained ellipsometric ratios plotted versus 1 / y/j gives a straight line from which 
the bending elasticity modulus can be determined. The success of the experiment not 
only demonstrates the correctness of the assumption of constant bending elasticity 
but also demonstrates the applicability of the capillary wave model. 

Exercises 

3.1 Derive eq.(3.1). 

3.2 Show that the capillary wave model in one dimension does not need the introduction 
of a bending elasticity modulus to obtain a mean squared thickness that is indepen- 
dent of a maximum wavenumber. 


9 see J. Van der Waals. Z. Phys. Chem. 13 (1894) 657 and J.W. Cahn and J.E. Hilliard. J. Chem. 
Phys. 28 (1958) 3219. 

10 see for instance B.P. Binks. H. Kellay. and J. Meunier. Thin Solid Films 210/211 (1992) 118-120. 
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4 Light scattering by single particles 

The scattering of light is described by van de Hulst 1 and by Kerker 2 . In these notes 
some important aspects of light scattering by single particles will be highlighted. 


4.1 Rayleigh theory of scattering by small dielectric spheres 


Consider a small dielectric sphere in the origin of a Cartesian coordinate system. The 
sphere has radius a and dielectric permittivity e 2 whereas the ambient medium has 
dielectric permittivity The dipole polarizability of such a sphere is 3 , see exercise 

4.1. 


a = 47ra 3 ei 


e 2 — e l 
62 + 2ei 


An incident electromagnetic field with electric field vector 


Ei — Eifie 


— ic ot 


(4.1) 

(4.2) 


induces in the sphere an oscillating dipole moment 


p = aE h0 e 


(4.3) 


Thus, the dielectric sphere acts as an electric dipole which radiates scattered waves in 
all directions. The wavenumber k = 2irni/\ where the refractive index n 1 = y/e[ and 
A is the wavelength of the light in vacuum. In a far point r, kr 1, the scattered 
electric field is given by 4 


gi kr 

E = k 2 (r x p) x r 

47reir 


(4.4) 


where r is a unit vector in the direction r , i.e. r = r/r. The direction of propagation 
is along r and the magnetic field vector is perpendicular to both the electric field and 
the direction of propagation. 

In what follows we shall assume that the incident light propagates along the z-axis 
of the coordinate system. The observer is located in the £Z-plane, which is then called 
the scattering plane , at a scattering angle 9 with respect to the z-axis, see fig. 4.1. 
For p-polarized incidence light, i.e. light that is polarized in the scattering plane, the 
scattered intensity is 


16t r 4 a 6 

p ~ 


e 2 — e l 
62 + 2ei 


2 

cos 2 9 


(4.5a) 


where L t = t\Ei ■ Ei* is the incident light intensity. For s-polarized light, i.e. light 
that is polarized perpendicular to the scattering plane, the scattered intensity is 

167r 4 a 6 


I, = h 


£2 - €1 


A 4 r 2 V.e 2 + 2ei 


(4.5b) 


^.C. van de Hulst. Light Scattering by Small Particles. Dover. NY 1957. 

2 M. Kerker. The Scattering of Light and Other Electromagnetic Radiation. Academic Press. NY 
1969, 

3 see: J.D. Jackson. Classical Electrodynamics. J. Wiley & Sons. New York 1975. sec. 4.4. 

4 see: J.D, Jackson. Classical Electrodynamics. J, Wiley & Sons. New York 1975. sec. 9.2, 
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4 LIGHT SCATTERING BY SINGLE PARTICLES 



Fig. 4.1. Scattering geometry: the observer is at r in the xz- plane (the scat- 
tering plane) and 9 is the scattering angle. 


Unpolarized incident light can be resolved into two incoherent linear polarized com- 
ponents which are parallel and perpendicular to the scattering plane. In this case 


Iu = 


Id + / s 


L 


LOTT a 


A 4 r 2 V.e 2 + 2e] 


(1 + cos 2 


(4.6) 


A random distribution of small identical dielectric spheres scatters light incoher- 
ently and the total scattered intensity is the sum of the contributions of the individual 
spheres: the scattered intensity is proportional to the number of scatterers. 

There are. in a historical perspective, two important consequences of the above 
given formulae: (1) the scattered intensity varies with 1/A 4 which accounts for the 
blue color of the sky and (2) for scattering angle n/2 the scattered light is polarized 
which accounts for the (partial) polarized sky light. 

In the case of absorbing spheres the dielectric permittivity is a complex number 
and the computation of the intensities changes likewise. 


4.2 Mie theory of scattering by spheres 

The general theory of light scattering by spheres is given by Mie in 1908. The solution 
can be written in terms of two amplitude functions, one for p-polarized light 


o _ V- 21 + 1 L dPl^ose) P}( cose) 

Sp{9) 2_ i / 1 , i \ i a i i/i T Ii 


l=1 l (l + 1) 

and one for s-polarized light 


d6 


sin 9 


2/ + i f P/(cosf)) dP}( cosS) 
bs{9) ~ \ a l 7T77, r Of 


/=1 l (l + 1 ) 


dd 


(4.7a) 


sin# 


(4.7b) 
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where P{ n {z) are the associated Legendre functions of the first kind of degree l and of 
order m 5 . The coefficients ai and bi can be conveniently expressed in terms of angles 
On an d Pi such that 


tan €q tan Pi 

di = : and bi = : 

tan ai — i tan Pi — l 

These angles, in turn, can be calculated by 


tana; 


nip[{y)pi{x) - ri 2 ipi(y)'ip'i(x) 

ni i’iiy)xi(x) - n 2 pi(y)xi(x) 


and 


n2p'i(y)Pi( x ) - n^ijyjipKx) 
" 2 ^i(y)xi(x) - n 1 pi(y)x'i(x) 


(4.8) 


(4.9a) 


(4.9b) 


with x = 2nn\a/ X and y = 2Ttn 2 a/X. The Ricatti-Bessel functions pi and \i are 
related to the spherical Bessel functions of the first kind ji and of the second kind, yi 

by 


Pi( z ) = z h{ z ) and Xi( z ) = z yi{ z )- (4.10) 

Great care has to be taken to evaluate the coefficients ai and 6/ numerically. It is for 
this reason that the coefficients are expressed in terms of the Ricatti-Bessel functions 
and not in the Hankel functions. Provided that one strictly follows the guidelines 
to evaluate these functions, as for instance put forward by Abramowitz and Stegun, 
the coefficients can be determined to any desired accuracy, see the appendix. The 
main reason to express the coefficients ai and bi in terms of the angles Q/ and Pi, cf. 
eq.(4.8), is to reduce the amount of complex number manipulations. 

For sufficiently large distances from the center of the sphere, i.e. kr 1 , the 
electric field of the scattered light for p-polarized incident light is also p-polarized 
with amplitude 

Ep = (4 - lla) 

and the electric field of the scattered light for s-polarized incident light is s-polarized 
with amplitude 

i -ifcr 

E s = -Ei—S 9 (e). (4. 

In fig. 4.2 some examples of scattering curves are given. The smallest latex sphere 
size (1 nm) is in the Rayleigh regime (the wavelength of the incidence light is cho- 
sen to be the HeNe line of 632. 8nm) and this illustrates nicely the cos 2 9 behaviour 
characteristic of p-wave scattering. The scattering for larger particle sizes tends to 
concentrate in the forward direction. 



5 M. Abramowitz and I.A. Stegun. Handbook of Mathematical Functions , Dover, NY 1964. 
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4 LIGHT SCATTERING BY SINGLE PARTICLES 



Fig. 4.2. Scattering intensities as a function of scattering angle for p-polarized 
incidence light of wavelength 632.8 nm. The particles are polystyrene latex 
spheres (n/=1.6) dispersed in water (n w —1.33). 


4.3 Rayleigh-Gans-Debije scattering 

Apart from the mentioned scattering functions for spheres there are also functions 
available of particles for other shapes such as stratified spheres and small ellipsoids 
and spheroids. These functions can be found in the monograph by Kerker. In this 
subsection we shall, rather then concentrate on a particular shape, discuss an ap- 
proximate theory for particles of arbitrary shape and size. The quantum mechanical 
analog of this theory is known as the Born approximation. 

The scattering function is calculated as follows: the particle, of arbitrary size 
and shape, is subdivided in infinitesimally small volume elements. Each such volume 
element is treated as a Rayleigh scatterer excited by the incident field which is sup- 
posed to be unperturbed by the prescence of the rest of the particle. For incident 
radiation polarized perpendicular to the scattering plane, the amplitude function for 
each volume element dV is given by 


dS s (6) 


. 3k 3 62 — 6i 

1 

47T e 2 + 2e l 


e i5 dV, 


(4.12) 


where exp(i<5) relates the phase of each elemental wavelet at the position of the 
observer to a common reference plane perpendicular to the incidence beam. At the 
observation point, far removed from the particle, the resultant amplitude function 
which arises from interference of each of the contributions is obtained by vector 
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summation 


s s (e) = i 


3 k 3 


e 2 - £i 


- J e iS dV. 


(4.13) 


47T e 2 + 2ei 

The fundamental approximation in this approach is that each volume element experi- 
ences the same incident field. For this reason neither the particle size nor the relative 
refractive index can become too large. For a relative refractive index n = sjc^/ei 
sufficiently close to unity eq. (4.13) simplifies to 


i k 3 


S s (9) = ^(n-1) J e iS dV. 


For real refractive index one then has for s-polarized light 

1 fP{9), 


k 4 V 2 

W) = ^(n 


47 r 2 r 2 

and for p-polarized light 
k A V 2 


h(») 


4i r 2 r 2 


(n — 1 ) 2 P(6) cos 2 (9). 


The quantity P(9 ) is known as the form factor and is given by 


m = fi\ f ^dV 


The form factor has been evaluated for a variety of configurations, 
himself gave the result for a homogenous sphere with radius a 

3(sinu — ucosu ) 2 


m 


9ir 


ir 


2m 3 J 3/2(«) 


(4.14) 

(4.15a) 

(4.15b) 

(4.16) 
Rayleigh 

(4.17) 


where u = 2kasm(6/2) and J 3 / 2 is the three-halves order Bessel function. 

Debije formulated in 1915 the theory for an aggregate of discrete elements rather 
than for a continuum 6 . If each of the scattering elements is identical and the aggre- 
gate assumes all random orientations, then the following average form factor can be 
calculated 


N N 


p(n\ _ * S1U d r ij 


(4.18) 


where is the distance between centers of the elements and q = 2k sm.(9/2). This 
expression has been evaluated for a number of configurations which are encountered 
in macromolecular chemistry. One of the most important applications being the 
randomly coiled high polymer chain for which Debije derived in 1947 

2(e~ w + w - 1) 


m = 


w * 


(4.19) 


where w = q 2 R \ and R g the radius of gyration of the polymer chain. 


6 see M. Kerker. The Scattering of Light and Other Electromagnetic Radiation. Academic Press. 
NY 1969. sec. 8.1, 
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4 LIGHT SCATTERING BY SINGLE PARTICLES 


4.4 Coupled-dipole method 

In the Rayleigh-Gans-Debije approach to the scattering of particles of arbitrary shape 
and size the main assumption was that the local field for each volume element was the 
incident field. This is no longer the case in the coupled- dipole method for calculating 
the light scattering by nonspherical particles 7 . The method consists of dividing the 
object into an array of N subunits located on a simple cubic lattice. The subunits 
are asumed to be sufficiently small to give rise to only electric dipole radiation. Each 
subunit is characterized by a (symmetric) polarizability tensor ay. such that the 
instantaneous dipole moment Pj can be expressed as 


Pj = a i\ E i-'52 A APk (4-20) 

\ k^j J 

where Ej is the electric field at the site j due to the incident light. The second 
contribution to the electric field A jkPk is due to the dipoles at positions k. 

Defining the matrix elements 

Ajj = aj 1 (4.21) 

one obtains the inhomogenous set of linear complex vector equations 

N 

Y, A *P*= E r ( 4 - 22 > 

fc=l 


The solution of this set of equations is a numerical problem and depending on the 
computational power at hand more or less success can be obtained. 

The important issue left to discuss is the actual number to put in for the dipole 
polarizabilities ay. The simplest approach is to use a Clausius-Mossotti or Lorentz- 
Lorenz relation (1.27) to calculate the polarizability. For spheres this does not always 
lead to satisfactory results compared to the exact Mie theory and one therefore often 
uses the electric dipole term of the Mie theory (4.8), i.e. 

6tt 

a = i-p-oq, (4.23) 

to obtain a better match. For nonspherical particles there is no way to assess the 
quality of the results apart from studying the convergence for increasing number N 
of subunits. 

An important problem of the method is the spurious enhancement of the scattering 
induced by the cubic lattice. Techniques have been proposed to damp out these effects 
by adapting the polarizabilities aj but. again, there is no way to assess the quality 
of the results for nonspherical particles. 


7 see C.E. Dungey and C.F. Bohren, J Opt Soc Am A 8 (1991) 81 and references therein. 
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Appendix: procedures to compute the scattering functions 

Below a set of Turbo-Pascal (Borland, Scott’s Valley, Ca.) procedures to calculate the Mie 
scattering functions for real refractive indices. 

• Definition of types, constants, and global variables. 


const 


l_max = 

10; 


lm_max = 

65; 


m_max = 

20; 

{ 


type 

real 

l_array 


2*l_max } 


= extended; 

= array [0..1_max] of real; 


lm_array = array [0 . . lm_max] of real ; 
var 


a_r ,a_i,b_r,b_i 

Q 

pie,tau 
Sl_r ,Sl_i 
S2_r ,S2_i 
11,12 


l_array ; 

real; 

l_array ; 

real; 

real; 

real; 


• Procedures to compute the Ricatti-Bessel functions 


procedure psi_calc(z : real; var psi,d_psi : l_array) ; 
const 

j_max = 100; 
var 

j : array [0..j_max] of real; 

j0,p : real; 

1 : integer; 

begin 

Miller’s algorithm utilizing recurrence relation backwards, 
see Abramowitz and Stegun, section 10.5 

> 

j [ j _max] : =0 ; 
j [j_max-l] :=1; 
for l:=j_max-2 downto 0 do 
j[l] : = (2*l+3)*j[l+l]/z - j [1+2] ; 
j0:=sin(z)/z; 
p:=j0/j [0] ; 
psi [0] :=j0*z; 
for 1 : =1 to l_max do 
begin 

j [1] :=j [l]*p; 
psi [1] :=z*j [1] ; 
d_psi [1] :=-l*j [1] +psi [1-1] ; 
end; 
end; 

procedure chi_calc(z : real; var chi,d_chi:l_array) ; 
var 

1 : integer; 

c,s : real; 

begin 

{ 

Simple forward recurrence relation algorithm, 
see Abramowitz and Stegun, section 10.5 

> 

c:=cos(z) ; 
s:=sin(z) ; 
chi [0] :=c; 
chi[l] :=c/z+s; 
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d_chi[l] :=c*(l-l/sqr(z))-s/z; 
for 1:=2 to l_max do 
begin 

chi [1] : =(2*1-1) *chi [1-1] /z-chi [1-2] ; 
d_chi[l] :=-l*chi [1] /z+chi [1-1] ; 
end; 
end; 

• Procedure to compute the functions tti(cos9) = (cos 6)/ sin 9 and ri(cos9) = 
do P^ (cos 9) 


procedure pie_tau__calc(z : real); 
var 

1 : integer; 

d_pie : l_array; 

h : real; 

begin 

{ 

Simple forward recurrence relation algorithm, 
based on recursion relation for associated 
Legendre functions with m=l. 

See Van de Hulst, sec. 9.31, Kerker sec. 3.7. 

Note: pi is predefined by Turbo Pascal 

> 

h:=l-z*z; 

pie[l] :=1; d_pie[l] : =0 ; 
pie [2] : =3*z ; d_pie [2] : =3 ; 
tau[l] :=z; 
tau[2] :=6*z*z-3; 
for 1:=3 to l_max do 
begin 

pie[l]: = ( (2*l-l)*z*pie [1-1] - l*pie[l-2] ) / (1—1) ; 
d_pie[l] : = (2*l-l)*pie [l-l]+d_pie [1-2] ; 
tau[l] :=z*pie[l]-h*d_pie[l] ; 
end; 
end; 

• Procedure to compute the coefficients a/ and bi. 


procedure ab_calc(ql ,q3 : real); 


var 

psil,d_psil,chil,d_chil : l_array; 

psi3,d_psi3 : l_array; 

tan_a,tan_b : real; 

1 : integer; 

h : real; 

begin 


Calculation of the coefficients a_l,b_l 
see Van de Hulst, p. 135 

> 

psi_calc(ql,psil,d_psil) ; 
psi_calc(q3,psi3,d_psi3) ; 
chi_calc(ql,chil,d_chil) ; 

Q : =0 ; 

for 1:=1 to l_max do 
begin 

tan_a:=-(ql*d_psi3 [1] *psil [1] -q3*psi3 [1] *d_psil [1] ) ; 

tan_a:=tan_a/(ql*d_psi3 [1] *chil [1] -q3*psi3 [1] *d_chil [1] ) ; 

tan_b:=-(q3*d_psi3 [1] *psil [1] -ql*psi3 [1] *d_psil [1] ) ; 

tan_b : =tan_b/ (q3*d_psi3 [1] *chil [1] -ql*psi3 [1] *d_chil [1] ) ; 

h:=tan_a*tan_a; 

a_r [1] :=h/(l+h); 

a_i [1] :=tan_a/(i+h) ; 
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h : =tan_b*tan_b ; 
b_r [1] :=h/(l+h); 
b_i [1] :=tan_b/(l+h) ; 

Q : =Q+ (2*1+1) *(a_r [1] +b_r [1] ) ; 
end; 

Q:=2*Q/(ql*ql) ; 
end; 

• Procedure that computes the scattering functions and intensities. 

procedure S_calc(ql,q3, theta: real); 
var 

1 : integer; 

z,f : real; 

begin 

{ 

Calculation of the amplitude functions 
see Van de Hulst, sec. 9.31 

> 

z:=cos(theta) ; 
pie_tau„calc(z) ; 
ab_calc(ql,q3) ; 

Sl_r :=0; Sl_i:=0; 

S2_r :=0; S2_i:=0; 
for l:=l_max downto 1 do 
begin 

f :=(2*1+1)/(1*(1+1)); 

Sl_r : =Sl_r+f * (a_r [1] *pie [1] +b_r [1] *tau [1] ) ; 

Sl_i : =Sl_i+f *(a_i [1] *pie [1] +b_i [1] *tau [1] ) ; 

S2_r : =S2_r+f * (b_r [1] *pie [1] +a_r [1] *tau [1] ) ; 

S2_i :=S2_i+f * (b_i [1] *pie [1] +a_i [1] *tau[l] ) ; 
end; 

II :=Sl_r*Sl_r+Sl_i*Sl_i; 

12 : =S2_r*S2_r+S2_i*S2_i; 
end; 


Exercises 

4.1 Consider an emulsion where equally sized droplets of the one fluid, dielectric constant 
eo, are suspended in another fluid of dielectric constant ei- Use the Clausius-Mossotti 
relation (1.22) and the expression (4.1) for the polarizability of the droplets to derive 
an expression for the dielectric permittivity e of a suspension with volume fraction p. 

4.2 Show that for small w the form factor (4.19) expands as 


P(9)e^ 0 = 1 - 


3 


(4.24) 


This leads to a possibility to measure the radius of gyration of a polymer by plotting 
P versus q 2 , the socalled Guinier plot . 



